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Abstract 

We use a configuration- interaction approach and Fermi golden rule to investigate electron-phonon 
interaction in realistic multi-electron quantum dots. Lifetimes are computed in the low-density, 
highly correlated regime. We report numerical evidence that electron-electron interaction generally 
leads to reduced decay rates of excited electronic states in weakly confined quantum dots, where 
carrier relaxation is dominated by the interaction with longitudinal acoustic phonons. 

PACS numbers: 73.21.La, 73.22.Lp, 73.61.Ey, 72.10.Di 



Semiconductor quantum dots (QDs) are artificially-confined systems where full control 
on the number of carriers and energy states can be achieved by either tuning the con- 

l|. These confined systems can reach unprece- 



finement potential or using external fields 
dented regimes of large Coulomb-to-kinetic energy ratios, t 



ius offering an ideal laboratory 



for the exploration of correlation effects at the nanoscale 3]. In addition, QD devices in 
the single-electron (SE) or multi-electron (ME) regime hold a primary role in new classes of 
applications, ranging from optoelectronics to biology and quantum information processing 

One fundamental question is how the discrete density of electronic states affects the re- 
laxation time of excited electrons. The knowledge of carrier lifetime is obviously of great 
importance for many applications, and most significant for the implementation of solid-state 
gates for quantum information processing, which require a large number of elementary op- 
erations to be realized within the phase coherence time of the relevant degrees of freedom. 
It was early recognized that, since scattering from optical phonon is suppressed in weakly 
bound QDs (with confinement energy below ~ 10 meV in GaAs) 0,0], the coupling with 
the longitudinal acoustic (LA) phonon continuum becomes the main source of scattering, 
thus limiting charge relaxation times to tens of nanoseconds |6|. Theoretical estimates of 
LA-phonon induced scattering rates were obtained in a SE approach and showed a satisfac- 
tory agreement with transport-spectroscopy experiments s7]. However, the technologically 
relevant problem of estimating the carrier lifetime when more than one electron are present 
in the device and, in general, the fundamental question of whether Coulomb interaction be- 



0. 



tween carriers affects the relaxation of ME quantum states have received little attention 

In this Letter we report numerical evidence that carrier-carrier Coulomb interaction leads 
to reduced relaxation rates of excited ME states in weakly confined, strongly correlated QDs 
when the dominating intra-band scattering mechanism is the interaction with LA-phonons. 
The lifetime for a specific QD is shown to strongly depend on the number of carriers. 

The calculation of ME states in QDs has been approached by several methods (for a review 
of quantum states in single and coupled QDs see Refs. QandO). S ince we are interested 
in the relaxation time of excited ME states, we need to know both ground and excited 
states with comparable accuracy. Our method of choice is the configuration interaction (CI) 
approach: the correlated wave functions of the ground and excited states are written as 
superpositions of Slater determinants (SDs), = rTa- c oj0)' obtained by filling in the SE 
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spin-orbitals a with the N electrons in all possible ways (c' a creates an electron in level a). 
The fully interacting Hamiltonian, written within the envelope function and effective mass 
approximations, is first block diagonalized, exploiting orbital and spin symmetries. Finally, 
the ME Hamiltonian is diagonalized in each symmetry sector, giving the full set of excited 
states of the system [loT ]. 

In order to evaluate quantitatively the relaxation time at zero temperature between given 
ME states b and a, n^ a , we employ the Fermi golden rule between fully correlated states jll|. 
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where the ME states \^ a ) = Y^i c ai\^i) have been written explicitly as linear combinations 
of SDs V q is the interaction operator of an electron with a LA-phonon with wave 

vector q, E a and hw q are the ME-state and phonon energies, respectively. We use the 
standard deformation-potential model for V q and a linear dispersion for the bulk LA-phonon 

HQ- 
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The SE levels are calculated for a cylindrical single or double QD, with confinement 
potential V(r) = V z (z) + ^m*^^ 2 + y 2 ), V z (z) representing the vertical confinement of 
a single or double quantum well, Huq being the SE energy level spacing of a lateral two- 
dimensional harmonic trap, and m* being the effective mass. We consider typical "vertical" 
GaAs/AlGaAs samples with effective cylindrical geometry For this configuration, 

and neglecting spin-orbit coupling, the SE spin-orbitals can be given the separable form 
i/>nmgcr( r , s ) — ¥>nm(%, y) ( t>g(. z )Xa ( s ) j with s =|, I the spin state, n = 0, 1, . . . the radial and 
m = 0, ±1, ±2, . . . the angular quantum numbers related to the lateral confinement, and 
g = 0, 1, . . . labeling the eigenstate in the z direction. Since the vertical confinement is 
much stronger than the lateral one, it is sufficient to consider g = for single QDs, and 
(7 = 0,1 for QDs coupled via vertical tunneling. In the ME Hilbert space, the resulting global 
symmetries are the total orbital angular momentum z-component M, total parity G, total 
spin S, and its z-component S z , which are used to block diagonalize the ME Hamiltonian. 
This numerical scheme has proved extremely reliable in reproducing the measured transport 
and inelastic light scattering spectra in single and coupled QDs |l0i \24\. 

Equation (0) shows explicitly that the ME scattering rate r^ a can be thought of as a 
combination of SE scattering events, where a given configuration (i.e., a single SD) changes 
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at most by the occupation of a single spin-orbital il> n mg<j( r , s )- The coefficients c a i are 
determined by the solution of the interacting problem through the CI method. It should 
be noted that the SE matrix elements ($j|V^|$j) depend strongly both on the SE orbitals 
ip nm (x,y)(j)g(z), possibly changing via confinement potential and external fields, and on the 
phonon wave vector q. The magnitude of q is determined by the transition energy E/, — E a , 
fixed, in turn, by the difference between the SE energies of initial and final levels plus the 
Coulomb energy. Indeed, typical LA-phonon wavelengths are comparable to the spatial 
extension of SE quantum states. Coulomb interaction, therefore, may affect r^_ ¥a through 
both the CI coefficients c a i and the SE matrix elements via shift of the transition energies 
Eb — E a . 

We start our discussion of ME-state relaxation times by considering the two-electron case. 
In Fig.QJwe show the relaxation rate r 6 l^ a (solid curve) for the LA-phonon induced transition 
as a function of the confinement energy ItjjJq in a single QD, where b = (M = 1, S = 0) is 
the correlated first excited state and a = (M = 0, S = 0) the ground state within the singlet 
sector. This transition is mostly composed of a SE p —>■ s transition [p = (m — 1, n = 0), 
s = {m = 0, n = 0)], with a "spectator" electron in level s (see diagram in Fig. [I}. For 
comparison, we show in Fig. ^ the corresponding non-interacting p — > s relaxation rate 
(dashed curve) |6j. The relaxation rate of the correlated state is suppressed by a factor 
~ 2.5 with respect to the SE transition. The two panels of the inset in Fig. ^ plot the 
same transition for the case of a double-dot sample with selected values of the inter-dot 
barrier. Here, the SE orbitals are coherently delocalized across the sample, with levels 
(mng) = (100) and (000) corresponding to p and s, respectively. The behavior is similar to 
the single QD system, with r 6 l >a vanishing at the low and high energy sides, but oscillations 
appear at intermediate values of fujo; these oscillations, which increase in number with the 



width of the inter-dot barrier, are a pure SE effect and are due to 



;he matching between 
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IHUJ. We stress 



the average separation of the two QDs and the phonon wavelength 
that, in the ME case, even if r b ~_^ a is reduced by electron correlation with respect to the 
corresponding SE transition, it behaves qualitatively in the same way. Specifically, the values 
of the confinement energy and of external fields at which the scattering rate is suppressed 
(cf. Fig. are the same. 

Now consider the lowest-energy SE transition in a parabolic QD, in which case the SE 
gaps hu>o fix the energy of emitted (bulk) phonons hw q = Huq. As the confinement energy 
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changes, e.g., by tuning the gate potential, both the phonon energy and the SE orbitals 
vary. In the SE case the scattering rate variation is due to the combination of these two 
effects. For the two-electron case, however, hw q is in general different from Hljq due to 
Coulomb contribution. Therefore, in order to verify whether r b ^ a is reduced by a genuine 
correlation effect or simply by the change in the emitted phonon energy, we calculated T^ a 
by artificially fixing fiw q at the SE gap fiuio also in the two-electron case. This is done by 
substituting the Dirac delta in Eq. (JTJ) with 5(fru — ftw q ). We found that the scattering 
rate obtained in this way (triangles in Fig. [TJ differs from the two-particle correlated case 
by only a few percents. Consequently, the reduction in r 6 l +0 is brought about by correlation 
only, since the transition energy and, therefore, the emitted phonon wave vector are only 
slightly affected by the interaction. 

We show next that the reduction of the scattering rate brought about by Coulomb cor- 
relation is a more general effect and holds also for a larger number of electrons. In the 
CI approach, the ME transition rate between two specific states is the square of a sum of 
many SE transition amplitudes [Eq. and, due to interference effects among those am- 
plitudes, there is no reason to believe a priori that the ME transition rates are larger or 
smaller than the SE ones. We found, however, that the calculated ME transition rate is, in 
general, smaller than the largest transition rate between two uncorrelated states of the QD, 
as discussed below. 

In order to discuss the scattering rates for N > 2, in Fig. El we show the LA-phonon 
induced scattering rate r 6 Z^ a as a function of the emitted phonon energy ftw q in QDs at 
selected lateral confinements and 2 < N < 6. Each point in the graphs represents a 
transition between two correlated ME states. The first 20 eigenstates are considered, for 
each value of N; due to total spin conservation, only a few of all possible transitions are 
allowed. In each panel several ME transitions accumulate precisely at hu>Q, identified by 
downward arrows: these transitions correspond to the de-excitation of Kohn modes Q, i.e., 
the separable center-of-mass modes of the ME system, whose energies are exact multiples 

Of ftjjjQ. 

To compare the ME scattering rates with SE ones, we report in Fig. El (dashed curve) the 
scattering rate of the lowest (p — *> s) SE transition for a QD whose parabolic confinement 
energy matches the emitted-phonon energy, as implied by the SE picture. Note that, for 
such a curve, the graph abscissa indicates also the confinement energy. Since each ME 
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transition has many SE components, in order to identify the transitions that are more 
directly comparable with the SE ones, we have indicated with the symbol X those processes 
whose highest-weight components of initial and final ME states can be mapped one into the 
other by destroying an electron occupying orbital p (initial SD) and moving it to s (final 
SD): the relaxation time of these specific ME transitions matches that of the SE p — > s 
relaxation once the correlation effects are neglected. 

At the low-energy side of each plot we see that a few ME transitions have a scattering 
rate which is higher than the SE one, which seems to contradict results for N = 2. However, 
note that here we are showing the ME scattering rates as a function of Ttw q at a given Fiui®, 
while in the SE case (dashed curve) the confinement energy is fojj qi as explained above. To 
single out the specific effects of Coulomb interaction on , a more significant comparison 
is with SE transitions (solid curve) calculated with the lateral confinement energy fuvo used 
for the ME calculations (values indicated in Fig. |2J), i.e. releasing the condition that equals 
the harmonic confinement to the phonon energy. In other words it is the solid curve of Fig. |2] 
that corresponds to the p —>■ s SE transition contributing to the ij summation of Eq. ^ As 
expected, the solid and dashed curves coincide at frw q = fkJo- For all considered cases, the 
scattering rate of ME transitions lies below the SE one with the same SE wave function and 
phonon energy; this is true, in particular, for Kohn modes. 

The three QDs considered in Fig. El correspond to rather different Coulomb-to-kinetic 
energy ratios in terms of A = lo/a%, where Iq = ^/H/m*uj is the QD radius and a% 
the effective Bohr radius, A = 3.3, 2.3, 1.9 for huo = 1,2,3 meV and GaAs parameters, 
respectively. The larger A (and the dot), the stronger the effect of correlation, keeping 
N fixed |20] . Figure El shows that larger QDs with stronger correlation effects have lower 
scattering rates. The above discussion demonstrates that Coulomb interaction tends to 
decrease the scattering rate of ME states in weakly confined QDs. From the point of view of 
specific applications, if a low SE LA-phonon scattering rate is predicted for a specific sample 
and/or field configuration, then such desirable property will not be spoiled by Coulomb 
interactions. 

The lifetime of a given quantum state depends both on the scattering rates and on 
the number of channels of possible de-excitation, i.e., the density of states of symmetry- 
allowed (spin-conserving) final states. In the ME system under study, an increasing number 
of particles leads to a larger number of states a given initial state can decay into. Such effect 
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competes with the above prediction of smaller scattering rates. We report in Fig. |3] inset 
the inverse lifetime of the first 20 QD excited states for fixed N, considering different curves 
for 2 < iV < 6. First we note that, due to the varying density of states, the energy range 
covered by the transitions considered decreases as the number of electrons increases. The 20 
decay times are reported, for each N, as a function of the difference between initial and final 
state energies, where the final state is the iV-electron ground state. The points are connected 
with lines as a guide to the eye. Although the curves do not show a monotonic behavior, 
two global features emerge: (i) as the energy of the initial state increases, its lifetime lowers 
(ii) the lifetime of a quantum state increases with N, provided the excitation energy is fixed. 
Property (i) holds since the density of de-excitation channels increases with the initial state 
energy, while (ii) is due to the fact that scattering rates decrease as correlation effects become 
stronger. In the samples we considered, the effect (ii) is not obscured by (i). This is shown 
in Fig. |3J where the scattering time averaged over the first 20 ME excited states is reported 
clS cL function of N for three different QD samples. 

To summarize, we quantitatively estimated the LA phonon-induced relaxation times of 
correlated ME excited states in weakly confined cylindrical GaAs/AlGaAs QDs. We find 
that electron-electron interaction contributes to reduce the scattering rate. Furthermore 
we showed that the correlated ME energy spectra mimic that of non-interacting QDs. In 
particular, the transition probability is suppressed for phonon energies whose wavelength 
matches the QD vertical dimension (when emission is inhibited also in the non-correlated 
case 14]), independently from the number of interacting electrons. We believe that these 
effects are of considerable relevance for the modeling of future nanodevices based on coherent 
electron dynamics and, in particular, for quantum computing proposals. 
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hco (meV) 



FIG. 1: Electron-LA phonon scattering rate at T=0 for: dashed curve, SE p — > s transition, 
and, solid curve, two-electron (M = 1, S = 0) — > (M = 0, S = 0) transition (see lower inset), 
for a cylindrical QD (GaAs layer 10 nm thick), as a function of the lateral confinement energy, 
fajjQ. Scattering rates for one and two electrons would be the same if correlation were switched 
off. Triangles show the two-electron rate calculated by artificially taking the energy of the emitted 
phonon equal to the confinement energy, as in the SE case: the values do not differ significantly 
from the "true" transition (solid line). The upper insets show the corresponding transition rates for 
a double dot with inter-dot barrier 3 nm (upper graph) and 8 nm (lower graph) wide, respectively. 
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FIG. 2: Electron-LA phonon scattering rate, for three QDs with different lateral confinement 
energy fajQ, as a function of the energy of the emitted phonon, hu> q . Dashed curves show the SE 
scattering rate of the lowest-energy transition (p — > s). Solid curves refer to the same transition 
except orbitals are fixed to the value of hu>o indicated in figure. Dots and X symbols identify specific 
ME transitions from excited to ground states with 2 < N < 6 (the first 20 excited states for each 
iV were considered), with the X symbols labeling ME transitions evolving into the non- interacting 
ones as correlation is switched off (see text). Arrows identify Kohn modes. 
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FIG. 3: Mean decay rate of correlated ME states with up to six electrons, due to electron-LA 
phonon scattering, for the three QDs of Fig.[2 with hw® = 1 meV (circles), 2 meV (squares), 3 meV 
(diamonds). Inset: decay rate of the first 20 correlated ME states as a function of the excitation 
energy for the QD with Hujq = 1 meV. The curves connecting the transitions considered are guides 
to the eye. The labels indicate the number of electrons. 
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